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We demonstrate that a simple phenomenological approach can be used to simulate electronic conduction in 
molecular wires under thermal effects induced by the surrounding environment. This “Landauer-Biittiker’s 
probe technique” can properly replicate different transport mechanisms: phase coherent nonresonant tunnel¬ 
ing, ballistic behavior, and hopping conduction, to provide results consistent with experiments. Specifically, 
our simulations with the probe method recover the following central characteristics of charge transfer in 
molecular wires: (i) The electrical conductance of short wires falls off exponentially with molecular length, 
a manifestation of the tunneling (superexchange) mechanism. Hopping dynamics overtakes superexchange 
in long wires demonstrating an ohmic-like behavior, (ii) In off-resonance situations, weak dephasing effects 
facilitate charge transfer. Under large dephasing the electrical conductance is suppressed, (iii) At high enough 
temperatures, UbT/cb > 1/25, with cb as the molecular-barrier height, the current is enhanced by a ther¬ 
mal activation (Arrhenius) factor. However, this enhancement takes place for both coherent and incoherent 
electrons and it does not readily indicate the underlying mechanism, (iv) At finite-bias, dephasing effects 
impede conduction in resonant situations. We further show that memory (non-Markovian) effects can be 
implemented within the Landauer-Biittiker’s probe technique to model the interaction of electrons with a 
structured environment. Finally, we examine experimental results of electron transfer in conjugated molec¬ 
ular wires and show that our computational approach can reasonably reproduce reported values to provide 
mechanistic information. 


I. INTRODUCTION 

Understanding of charge transport mechanisms in 
single-molecule junctions is essential for the realization 
of molecular electronic devices, as well as for eluci¬ 
dating many processes in chemistry, biology, and con¬ 
densed phase physics. Examples include charge transfer 
in DNA 1,2 , a process which plays a crucial role in muta¬ 
genesis and carcinogenesis, electron transfer reactions 3-5 , 
and electron correlation effects at the nanoscale, such as 
in quantum dots and molecular junctions 6-8 . 

Experiments probing electron transfer (ET) rates 
in donor-acceptor DNA molecules 9-12 , self-assembled 
monolayers and single-molecule wires 13-31 have demon¬ 
strated the central role of two limiting transport mech¬ 
anisms: phase coherent nonresonant tunneling (superex¬ 
change) and incoherent thermally-activated hopping. Di¬ 
rect tunneling, a single step transition, dominates trans¬ 
port in short wires, but becomes unlikely at large molecu¬ 
lar lengths as tunneling rates decrease exponentially with 
length. In contrast, long-range electron transfer can take 
place when the transport process is broken into multiple 
steps, with the electron (or hole), transiently localized on 
molecular sites, hopping between them. Furthermore, in 
molecular wires bridging two voltage-biased electrodes, 
“field emission”, (Fowler-Nordheim tunneling) behavior 
develops when the applied voltage bias transforms the 
tunneling barrier from a rectangular form into a triangu¬ 
lar shape 20,24 . 
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In an effort to explain experimental results, and 
moreover, predict molecular electronic functionality, a 
plethora of theoretical and computational methodologies 
were developed, aiming to explore the role of environmen¬ 
tal effects (internal molecular motion or the surround¬ 
ing matrix) on molecular conduction. These tools can 
be roughly grouped into two classes: (i) microscopic- 
physical modelings which are valuable for small systems, 
and (ii) phenomenological descriptions, compromising 
the completeness and exactness of the model to enable 
large-scale calculations. The first type of approaches cor¬ 
responds to models in which molecular vibrations, and 
other many-body interactions such as electron-electron 
repulsion are explicitly included in the model Hamil¬ 
tonian. The dynamics and steady-state properties of 
the system can then be analyzed by a variety of treat¬ 
ments. A non-exhaustive list, focusing on vibrational 
effects in molecular conduction, includes density oper¬ 
ator approaches 32,33 , Green’s function tools 34 and path 
integral simulations 35-37 . Since the vibrational degrees 
of freedom and/or electron-electron interactions are ex¬ 
plicitly incorporated in such treatments, simulations are 
restricted to minimal models with a single molecular elec¬ 
tronic state, or a pair of states. The Anderson-Holstein 
model takes into account a single electronic level and a 
particular vibration. Other common modelings consider 
only two molecular orbitals, HOMO and LUMO, as de¬ 
cisive for electronic conduction 34 . 

Complementing minimal-microscopic approaches, it 
is beneficial to establish effective-phenomenological 
methodologies for electronic conduction with a more fa¬ 
vorable scaling with molecular size. The Lindblad for- 
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FIG. 1. (a) A molecular wire with N electronic sites (inter¬ 

site tunneling v and bridge height es) coupled to a thermal 
environment, represented by the shaded region surrounding 
the electronic states. The thermal bath, including molecular 
and external vibrational degrees of freedom, is characterized 
by its temperature T, spectral density function with a cutoff 
frequency fl and electron-vibration coupling strength 7^. (b) 
In this work we introduce the environment in a phenomeno¬ 
logical manner by using the probe technique, by attaching 
fictitious metal terminals to electronic sites. The figure illus¬ 
trates voltage probes, characterized by their temperature and 
chemical potentials which are determined such that there is 
zero net charge current to each probe. 

malism, Redfield theory, and other kinetic equations, 
are examples of such approaches 5,38 . In these treat¬ 
ments the interaction of electrons with environmental 
degrees of freedom is incorporated into the dynamics 
via decoherence and dissipation rates which are intro¬ 
duced into the equations of motion for the reduced den¬ 
sity matrix. Variants of such projection operator tech¬ 
niques have been employed extensively in the literature 
to model electron, proton, phonon, and exciton dynamics 
in condensed phases 4,5 . In donor-bridge-acceptor elec¬ 
tron transfer (DBA-ET) processes, such methods have 
proven useful for describing, in a unified manner, both 
the tunneling and the incoherent-hopping regime 39-48 . 
More recent works rigorously examined the consistency 
of these effective treatments for describing steady-state 
characteristics 49-51 , and applied kinetic approaches to 
predict electronic functionality 52,53 . 

A distinct phenomenological route to implement deco¬ 
herence and inelastic effects was introduced by Biittiker 
in Ref. 54 . In this “Landauer-Biittiker’s probe” (LBP) ap¬ 
proach the non-interacting Hamiltonian is augmented by 
probe terminals in which electrons lose their phase mem¬ 
ory and exchange energy with probes’ degrees of free¬ 
dom. This technique, so far mostly employed in meso¬ 
scopic physics 55 , e.g. to study heat to work conver¬ 
sion efficiency in thermoelectric devices 56-58 and mag¬ 
netotransport in quantum dot systems 59,60 , is particu¬ 
larly appealing: It allows one to model phase breaking 


processes e.g. in quantum dots, while using the Lan- 
dauer non-equilibrium Green’s function (NEGF) com¬ 
putational tool 6-8,61 . Beyond electronic conduction, the 
LBP technique has been adopted to explore heat trans¬ 
port and thermal rectification in linear chains 62-68 and 
two-dimensional constrictions 69 . 

The LBP method had proved itself predictive in studies 
of charge and heat transport in mesoscopic physics. Can 
it appropriately describe environmental effects on elec¬ 
tronic conduction in molecular junctions? Recently, this 
technique has been applied to explore the role of decoher¬ 
ence on quantum interference in molecular conduction' 0 
and to examine the connection between charge trans¬ 
fer kinetics and steady state currents 71 , showing that 
it can provide both the tunneling mechanism for short 
bridges, and a “soft”-ohmic distance dependence for 
longer chains. However, a careful exploration of the LBP 
method in molecular electronic configurations is missing. 
Our objective here is to employ the LBP technique and 
examine whether it can yield results consistent with our 
comprehensive knowledge of environmental-assisted elec¬ 
tron transfer effects in molecular electronic junctions. 
Particularly, the technique should depict the turnover 
in transport, from tunneling to hopping, with increasing 
molecular size, temperature, and dephasing rates. 

Another objective of this work has been to develop 
an approximate analytic expression for hopping conduc¬ 
tion in molecular junctions. Several works have un¬ 
dertaken this task in a donor-bridge-acceptor configura¬ 
tions to construct the electron transfer rate , a related 
quantity 42 ’ 44-4 '’' 2 ’' 3 . However, few theoretical studies 
had considered this problem in a metal-molecule-metal 
setup 28 ’ 43,74-76 , which differs fundamentally from the 
DBA-ET situation: While in the latter case the initial 
condition places all electrons in the donor state with 
a given energy, in a metal-molecule-metal experiment 
many electronic states in the metals contribute, and elec¬ 
trons, with a thermal distribution of energies, determine 
the electronic conductance. 

Simulations presented in this work confirm that 
the LBP method can properly emulate experimentally- 
observed characteristics of molecular electronic conduc¬ 
tion in wires with a tunable “dephasing strength” param¬ 
eter, characterizing the strength of electron-environment 
interactions. Based on our simulations we suggest an 
analytic expression which describes hopping conduction 
in molecular wires. The Landauer-NEGF technique is 
nowadays the leading method in modeling the opera¬ 
tion of phase-coherent molecular electronic devices 6-8,61 . 
Particularly, recent measurements of molecular ther¬ 
mopower (see for example Refs. 7 '’ 78 ) and heat dissipation 
in nanojunctions 79 were explained using this method. By 
establishing here the appropriateness of the LBP treat¬ 
ment in molecular electronic conduction problems we ad¬ 
vocate for the natural generalization of the Landauer- 
NEGF technique to incorporate environmental-assisted 
transport effects by including LBPs. 

This paper is organized as follows. In Sec. II we 
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present the wire model and the probe technique. Sec. 
Ill summarizes transport mechanisms in molecular junc¬ 
tions. In Sec. IV we present numerical simulations under 
low applied voltage. Based on observations, we construct 
an analytic expression for bath-assisted molecular elec¬ 
tronic conduction. The behavior of the junction under 
large applied voltage is described in Sec. V. In Sec. VI we 
apply our computational method and compare findings 
to experimental results for transport in long conjugated 
molecules. We conclude in Sec. VII. 

In this work we interchangeably refer to the molecular 
system as a “molecular wire”, “molecular junction”, or, 
a “molecular bridge”. Identifying the system as a “wire” 
emphasizes that molecules under consideration include 
repeating units, and that they are capable of transfer¬ 
ring charges over long distances. The term “junction” 
highlights the geometry: the molecule is placed between 
two metal electrodes, distinguishing it from the DBA sit¬ 
uation. The title “bridge” describes the structure as¬ 
sumed: We study situations in which molecular orbitals 
relevant for transport lie above the Fermi energies of the 
two electrodes. Molecular electronic degrees of freedom 
are coupled to environmental coordinates (here, mim¬ 
icked by probes). The environment may correspond to 
internal molecular vibrations or the solvent’s degrees of 
freedom (we exclude the metals’ electrons from this defi¬ 
nition). We collectively refer to these nuclear coordinates 
as a “thermal bath”, “environment”, or “surrounding”, 
assumed to be maintained at the temperature of the met¬ 
als. 


II. MODEL AND METHOD 
A. Model Hamiltonian 

We consider a molecule bridging metal electrodes with 
the Hamiltonian 

H = H w + H l + H r + H t + H p + V p . (1) 

It comprises the molecular wire Hw , two electrodes Hl 
and Hr 1 and a coupling Hamiltonian Hp allowing charge 
transfer between the two leads and the wire. We further 
introduce (in a phenomenological manner) dephasing and 
inelastic effects for electrons on the molecule. This is 
achieved by attaching local reservoirs (Hp) to each site 
in the wire (coupling Hamiltonian Vp). see Fig. 1. The 
molecular wire includes N single-level sites of energy e n , 

N N -1 

H-W = ^ ' CncljCn T ^ ( Un,n-j-lC^Cn-l-l “b h.C. (2) 

n =1 n —1 

Here, c 4 (c n ) are fermionic creation (annihilation) oper¬ 
ators of electrons on each site in the wire, the parame¬ 
ters v n , n+ i are the inter-site tunneling energies. In what 
follows we consider molecules made of identical building 
blocks, thus we introduce the short notation ep = c n and 


v = Vn,n+i • The metal electrodes are modeled by a Fermi 
sea of noninteracting electrons, 

= y ^ = T, R. (3) 

k 

k (a K k) are fermionic creation (annihilation) operators 
of electrons with momentum fc in the v lead. Electrons 
can tunnel from the L ( R ) metal to site 1 (V), 

Hp = E 9L,ka{ k C! + ^2 9R,k^ Rk c N + h.c. ( 4 ) 

k k 

In the absence of the probes, this Hamiltonian dictates 
phase-coherent electron dynamics, reflected e.g., by a 
tunneling behavior 5 . We now include N probes, addi¬ 
tional metal electrodes 

N 

Hp = 'y ) y ) £n,/cn]j (5) 

n— 1 k 

The nth probe can exchange particles with the nth site 
of the molecular wire, 

N 

Vp = E E 9 n,kal k c n + h.c. (6) 

n=1 fc 

Here a' k ( a nt k ) are fermionic creation (annihilation) op¬ 
erators for an electron in the n = 1, 2,..., N probe with 
momentum fc, g n ^ are the tunneling energies from the 
nth molecular site into the nth probe. To eliminate 
charge leakage processes from the molecular wire into the 
probes we enforce certain conditions on conduction. In 
Sec. IIB we describe in detail two such constrains, the 
“dephasing probe” and the “voltage probe”. 


B. Biittiker’s Probes technique 

The Landauer approach provides a simple-exact de¬ 
scription of phase coherent quantum transport 6 . Given 
its simplicity, it is appealing to use it beyond the coherent 
limit. Indeed, as was shown in Ref. 54 , one can implement 
(elastic and inelastic) scattering of electrons with other 
degrees of freedom, possibly phonons, photons, and other 
electrons, by introducing additional terminals (probes) 
into the model system. The key point here is that the pa¬ 
rameters of these terminals, essentially their local charge 
distributions, should be set in a self-consistent way such 
that there is no net (average) particle current between 
the physical system of interest and the probes. 

The probe technique can be implemented under differ¬ 
ent self-consistent conditions, allowing us to craft elec¬ 
tron scattering processes: elastic dephasing effects are 
implemented via the “dephasing probe”, while dissipa¬ 
tive inelastic effects are introduced through the “voltage 
probe”. Further, dissipation-less inelastic scattering pro¬ 
cesses can be admitted by requiring the net average fluxes 
of particles and heat between the probe terminal and the 
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system to vanish, termed as the “temperature-voltage 
probe”. These probes can be operated in the linear re¬ 
sponse regime, as well as far from equilibrium 80 . 

Below we employ v to identify the L and R (physical) 
metal electrodes to which the molecule is connected. We 
count the probe terminals with the index n and use a to 
identify all leads, the two metal electrodes v = L,R and 
the n = 1,2,.., TV probes. 

Since the model Hamiltonian does not include interac¬ 
tions, its charge transfer characteristics can be described 
with the Landauer-Buttiker formalism 6 . The total cur¬ 
rent leaving the L contact is given by 

_ f-OO 

Il = ^h Y, J oo 7i, a (e)[/ L (e)-/ a (e)]de. (7) 

This expression should be multiplied by a factor of 2 to 
account for the spin degree of freedom. Since magnetic 
effects are absent, T a , a '{e) = 7^', a (e). /„(e) = [e^^~^G-\- 
1 ] _1 are the Fermi functions in the physical electrodes, 
given in terms of the inverse temperature k R T = 
and chemical potentials The functions f n (e) are to 
be determined from the probe condition. 

The electrical conductance is defined as the ratio of 
charge current to applied voltage, AH = (/zl — g i#)/e, 

G = I l /AV. (8) 

This definition does not necessarily assume low voltages. 
One often further defines a linear response electrical con¬ 
ductance demanding the applied bias to be the smallest 
energy scale in the system, |eAH| < 1//3, D, 7 „, es, v, 7 d, 
with 7 j, and 7 ^ defined after Eq. (13) and D as the band¬ 
width of the metal leads, the largest energy scale. In our 
simulations we principally work under the low-bias con¬ 
dition, yet for simplicity use everywhere the definition 
( 8 ). 

In direct analogy to Eq. (7), the net current between 
the nth probe and the system can be written as 

_ pOO 

In = J Tn ’ a ^ de ' ^ 

The transmission functions in Eqs. (7) and (9) are ob¬ 
tained from the (TV x TV) Green’s function and the hy¬ 
bridization matrices 5 , 

T a , a .(e) = Tr[r a (e)G^(e)r ct /(e)G(e)], (10) 

where the trace is performed over the TV states of the 
molecule. The Green’s function is obtained from the in¬ 
verse of a tridiagonal matrix, A = [Gi] -1 , with the ma¬ 
trix elements 

A„. n (e) = e - e B + -[7l(<0<W + lR^)^n,N + 7™(e)] 

A„ >n ±i(e) = -v. ( 11 ) 

The hybridization matrices have a single nonzero value 

as follows, 

[I7n(c)]n,n — 7ro(^)> 

[r L (e)]i,i= 7 i(e), [r_R(e)]AT,Ar = 7 /i(e), ( 12 ) 


with energies 

7<*(*) = 2 tt ^ \9a>,k\ 2 8(e - €«,*). ( 13 ) 

k 

We work in the wide-band limit, unless otherwise stated, 
and take as energy independent parameters. Below 
we assume that all sites in the wire are similarly affected 
by the probes, thus we use a single parameter to identify 
the probe-molecule hybridization energy, 7 d = r y n - We 
also define the harmonic mean 7 = ^7_l7.r/(7l + 7_r) as 
a measure for the molecule-metal coupling. 

We now describe the dephasing probe and the voltage 
probe conditions. In the linear response regime at low 
enough temperatures the two probes act similarly on the 
molecule. 

Dephasing Probe. We implement local elastic dephas¬ 
ing effects by forcing the energy-resolved, particle cur¬ 
rents, between every probe and the system, to diminish. 
Mathematically, we demand the integrand in Eq. (9) to 
nullify. This condition translates into a set of TV linear 
equations for the functions / n (e), given in terms of the 
Fermi distributions at the L and R electrodes, 

+ r -T(e) + T n , R {e) fn(e) ~ ]T T n , n >{e)f n >{e) 

_ n' J n' 

= T n , L {e)f L {e) + T n , R {e)f R (e). (14) 

We solve this linear set repeatedly, for energies e extend¬ 
ing —D to D , with D a high energy cutoff (bandwidth). 
The resulting probe distributions /„(e), not necessarily 
in the form of Fermi functions, are used in Eq. (7) to 
reach the net charge current. 

Voltage Probe. Local inelastic-dissipative effects can 
be introduced into the wire by demanding that the net 
total particle current flowing between each probe and the 
system, Eq. (9), vanishes, 

In = 0 . ( 15 ) 

Out-of-equilibrium, these TV equations can be solved nu¬ 
merically, as was done in phononic models 68 . Analytic 
results for p n can be reached in linear response 56 , 

^rs~ a -th Mik 

= J de [T n ,L(e)lJ.L +T n ,R(e)fj, R ]. ( 16 ) 

We can further simplify Eq. (16) at low enough temper¬ 
atures, 1//3 < 7j/. In this limit, the transmission function 
is assumed to be flat in the vicinity of the Fermi energy, 
in a region where the Fermi function suffers significant 
changes, 21^1 ^ —S(e — e_p). 

At low applied bias, the dephasing and voltage simi¬ 
larly operate. This can be justified by noting that in lin¬ 
ear response and at low enough temperatures the voltage 
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probe can support energy exchange processes only within 
a small energy interval, /i n — /if, thus it acts similarly to 
the dephasing probe (see Fig. 8 ). 

What does a dephasing probe do? While in the local- 
site basis this probe enforces the condition of zero en¬ 
ergy resolved flux between the every site and its lo¬ 
cal probe, upon transforming the Hamiltonian of the 
wire ( 2 ) from the site representation to the energy ba¬ 
sis we find that the dephasing probe can drive transi¬ 
tions between molecular orbitals: We define the molec¬ 
ular eigenstates by the fermionic creation operator 
Hw = Em Emdlndm, then formally expand the local- 
site operators, c n = E m V™^m- In the energy basis 
the N probes are coupled to the N molecular orbitals, 
Vp = J2n,m^n,m9n,kdl l k dm +h.c ., providing a mecha¬ 
nism for (environmentally-mediated) transitions between 
molecular electronic states. The energy resolved current 
to the n probe is given by 

In,k = ' ( q 9n,k^n,m ( ^rnb'n,k — 9n,k^n,md\ l k d rr ^j) 

m 

(17) 

Since under the dephasing probe we demand that = 
0 , this expression reveals that the n probe takes electrons 
of energy ek e.g. from the mth orbital and disperses them 
into the other orbitals. For example, If we consider a 
single probe and two molecular orbitals, mi and m 2 , the 
dephasing probe condition translates to the conservation 
law I rni —tprobe (e) = I pro be^m 2 (e ), he., the net energy flow 
to the probe is zero, but it exchanges electrons between 
orbitals. 

Below we refer to the hybridization energy 7 ^ 
(molecule-dephasing probes) as “dephasing strength”. 
As explained above, 7 d/h does not correspond to a “pure 
dephasing” process. This situation is similar to that 
reached with the perturbative (system-bath weak cou¬ 
pling) Redfield equation 43 : Bath-assisted dynamics had 
been admitted to the model by adding dephasing oper¬ 
ators, acting locally on different sites of the wire. In 
the energy basis of the molecule these local bath oper¬ 
ators were responsible for bath-induced transitions be¬ 
tween molecular orbitals 42 . Given the probe conditions, 
the current leaving the L contact is identical to the cur¬ 
rent reaching the R terminal, Ip = —Ir. 


C. Structured environment 

The interaction of electrons with molecular vibrations 
is explicitly described by the following model 


H = H w + H L + H R + H T + H ph + He- ph . (18) 

The molecular wire and the electrodes are the same as in 
Eqs. (2)-(4). The bosonic (phonon) environment and its 


coupling to electrons is given by 


Hph — 


^ ' bJq,nb^ q n bq t n, 
q,n 


H P - 


ph 


y ' hq, n (b q n + bq ^C^Cn- 

q,n 


Here fef n (b qn ) stands for a bosonic creation (annihi¬ 
lation) operator responsible for exciting a mode of fre¬ 
quency q , coupled to the electron number operator c' n c. n . 
The environment and its coupling energies h qtU provide 
the “spectral density” function J„(w) = |/i gjrl | 2 c>(w — 
w q , n ), possibly comprising a broad (low frequency) con¬ 
tribution, as well as discrete high frequency vibrational 
modes. It was illustrated in many studies that structured 
environments support involved dynamics in the subsys¬ 
tem; essentially, the back action of the bath on the system 
drives “non-Markovian” dynamics 38 . 

To model such effects within the LBP technique 
we suggest to adopt energy-dependent probe relaxation 
rates, defined in Eq. (13). Specifically, we use below the 
Debye-Drude function 

7d(e) = 27d e2 ^ 2 (19) 

which reaches the maximum value 77 at the cutoff fre¬ 
quency, 7 d(e = fi) = 7 d- Electrons arriving at the nth 
electronic site with energy |e| < H are likely to exchange 
energy with the probe’s degrees of freedom. In con¬ 
trast, electrons will not scatter to the probe terminals if 
|e| 3> fl, and they will coherently cross the junction. This 
behavior corresponds to the expected effects of electron- 
phonon interactions. In Sec. IV F we examine the role of 
a finite fl on the conductance. 


III. TRANSPORT MECHANISMS 

In this Section we summarize results for the electrical 
conductance of molecular wires in (i) the coherent limit 
when 7 d = 0, see Fig. 2(a), and (ii) assisted by the 
environment with 7 ^ ^ 0, see Fig. 2(b). 


A. Coherent electrons: Ballistic motion and Tunneling 

In the absence of the probes, 7 d = 0, the current obeys 
the standard Landauer expression, 

II = ^ deT^nmUe) - / fl (e)]. (20) 

Pronounced temperature dependence of transport char¬ 
acteristics is typically attributed to the hopping mecha¬ 
nism, but the current can also strongly depend on tem¬ 
perature in coherent tunneling due to the Fermi functions 
in this expression. We illustrate this next. We consider 
the contribution of a single-dominant narrow resonance, 

TL,fl(e) ~ 'yA(e B )6(e - e B )- (21) 
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FIG. 2. Dominant transport mechanisms in the LBP molecular wire model, (a) Off-resonant tunneling and (resonant) ballistic 
conduction when 7 d = 0. The blue line at the left end describes the Fermi distribution function of incoming electrons, (b) 
Illustration of bath (probe) assisted transport, depicting for simplicity a single probe. While all molecular orbitals are broadened, 
we explicitly demonstrate the broadening of two orbitals, those that are strongly coupled to the left and right metals. Incoming 
electrons of energy e, possibly at the tail of the broadening function, are scattered between the probes and the molecular states, 
then emitted to the R lead. Energy conservation is imposed by the dephasing probe. 


Here A as a dimensionless prefactor. We highlight its 
dependence on the bridge energetics, but it may further 
depend on other molecular parameters, inter-site tunnel¬ 
ing v, hybridization to the electrodes 7^ and molecular 
length N. The prefactor 7, the harmonic average of 
and 7#, has been included here as the resonance width. 
Using this narrow-resonance form, Eq. (20) reduces to 

II = ( e s)[/i( e s) - /it(es)] 

= A(e B )e^[e-^-e-^]f L (e B )f R (e B ). 

( 22 ) 

In the linear response limit, /3A/i < 1, we expand the 
functions to the lowest nontrivial order in voltage, 

G = Go7/U4(e B ) + 

G 0 p'yA(e B )e-^ B . (23) 

The combination /3e suppresses the tunneling con¬ 
ductance at both low and high temperatures, provid¬ 
ing a temperature dependence distinctively different from 
standard Arrhenius behavior. Here Go = e 2 /27 tH; we left 
out the factor of x 2 coming from the spin degree of free¬ 
dom. 

Going back to Eq. (22), we now consider a finite bias 
situation, /3A/i > 1. In this case we maintain only the 
larger factor in the numerator of Eq. (22), for A/r > 
0, and again take the limit /3e B 3> 1- We now reach an 
Arrhenius activation form for the conductance 26 , 

G = Go 7 A(e B )e _/3(ei3 " |eAy|/2) . (24) 

Eqs. (23)-(24) describe contribution to conductance from 
thermally excited electrons, populating energies in res¬ 
onance with a narrow molecular orbital of energy e B . 
These electrons cross the molecule in a ballistic manner 


as dephasing and inelastic terms are missing from these 
formulae. 

Next, we broaden the resonance (21). This allows 
electrons at the Fermi energy to cross the bridge, off- 
resonance. We assume (i) that the applied bias is small 
to expand f v (e) ~ / eg (e) - ~ t F ), and (ii) that 

the transmission function is almost flat (constant) close 
to the Fermi energy. We then reach the following, well- 
known, result 

lL = A ^ TL ’ R{eF) * 

= ^t^AA^f)- ( 25 ) 

The conductance is proportional to the transmission eval¬ 
uated at the Fermi energy, and it does not depend on 
temperature, 

G = GoTlA^f)- (26) 

In the “deep tunneling” regime when the barrier height 
is large, e B > the transmission function can be 

simplified into the “superexchange” expression 5 

7L,fl(eF) ~ 

IlIrv 1 

[A e F ) 2 + ( 7 l/2 ) 2 }[(e B - e F )* + ( 7 b /2 ) 2 ] 

/ \ 2 (N— 2 ) 


\e B — e F J 

/ \ 2N 
cb » 7 *, 7 LlR ( v \ 

v 2 \e B -e F ) 



Recall that we set energies relative to the Fermi energy 
(e F = 0). Eq. (27) indicates that the conductance decays 
exponentially with molecular size, L = aN , with a as a 
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unit length. Identifying the tunneling decay constant by 
k = | In | ^ | (often denoted by f3 in the literature, here 
reserved for the inverse temperature), we arrive at the 
familiar form 


G ~ Goe~ KL . (28) 

Fig. 2(a) depicts the ballistic [Eqs. (23)-(24)] and tun¬ 
neling [Eq. (27)] contributions to the current. 


B. Bath-assisted transport 

Calculations based on perturbative master equation 
formalisms suggest that long-range electrical conduction 
may be captured by a rational function 43 , correspond¬ 
ing to experimental observations of long-range ohmic-likc 
conduction in molecular wires 9,10,12-15,18-22 . In Sec. IV 
we present detailed simulations using the LBP technique, 
and based on observations, we construct the following 
form for the hopping (H) conduction, 

Gk-WT)*-^. (29) 

It is valid when the bridge is high es f, 7 d, A/x, as well 
as e_e > 7, and in the range 15 < f3es < 30. This cor¬ 
responds to cb = 0.5 eV and temperatures in the range 
±150 K around room temperature. As before, N stands 
for the number of repeating units in the wire, 7^ is the de¬ 
phasing strength, v the intersite tunneling, es the bridge 
energy relative to the Fermi energy, A(T) a dimension¬ 
less prefactor which depends on temperature, possibly 
in a form weaker than the “bare” Arrhenius factor, e.g., 
A(T) ~ /3e~ l3EA , with an activation energy Ea , linearly 
proportional to cb■ The parameter l is introduced to 
accommodate the tunneling-to-hopping transition point, 
see for example Fig. 3. 

Similarly to environmentally assisted DBA-ET rates, 
Eq. (29) describes an ohmic conduction, it depends on 
v 2 , and as expected, grows with 7 d- However, an ex¬ 
act correspondence to ET rates is missing 41,42,44-47,72,73 . 
While these two quantities, DBA-ET rate and the elec¬ 
trical conductance (examined here) are obviously related, 
their properties are fundamentally different in some ways: 
DNA-ET processes involve electron localization on the 
donor and acceptor sites whereas in a molecular wire 
setup electrons are transferred through the bridge be¬ 
tween metal electrodes. Thus, the hybridization energy 
of the bridge to the metals and the thermal occupation 
factor of electrons in the bulk largely determine the elec¬ 
trical conductance. Furthermore, in DBA situations the 
Franck-Condon factor, taking into account the energy 
difference between donor and acceptor states and the 
solvent-induced reorganization energy, dictates the trans¬ 
mission rate, while it is missing from molecular electronic 
conduction 75,76 . Theoretical studies have predicted a lin¬ 
ear relationship between DBA-ET rates and the conduc¬ 
tance (for molecules with high barriers), but experiments 


indicate that these two quantities are not simply linearly 
correlated 81 , possibly due to differences in bath-induced 
decoherence rates and bridge energetics 71 . 

Fig. 2(b) schematically depicts the probe-assisted hop¬ 
ping conduction, Eq. (29). While in DBA-ET situa¬ 
tions electrons can populate the bridge only if thermally 
excited, in the junction setup the molecular orbitals 
are broadened due to their hybridization to the metals 
[Lorentzian lineshapes in Fig 2(b)], and this broadening 
allows electron to occupy molecular states without ther¬ 
mal activation. 

We further argue that at low temperatures, equation 
(29) should be generalized to 

N / rn—\ \ ^ 

G D = G 0 l^r) F m{lL,lR,ld,T). (30) 

m=2,3.. \ e B J 

interpolating between the coherent result (27) using 
Fn = IlIr and F m <jv = 0 and the hopping limit (29) 
when F 2 ~ A(T)j 2 /(N + l) » F m>2 . 



FIG. 3. Electrical conductance as a function of molecular 
length N at T = 300 K using different dephasing strengths, 
7d=0, 1, 10 and 100 meV, as indicated in the figure cap¬ 
tion, (dephasing probe condition), (a) A semi-logarithmic 
plot demonstrates that G decays exponentially with length in 
short chains at weak dephasing, (b)-(e) An inverse algebraic 
behavior (29), G -1 oc N, is established for long chains at 
large enough dephasing. Dominating transport mechanisms 
are marked in panel (a). 


IV. RESULTS: TUNNELING TO HOPPING 
CONDUCTION 

We describe here the LBP numerical simulations of 
molecular electronic conduction (mostly under the de¬ 
phasing probe), with the goal to construct an approx¬ 
imate closed-form expression for the hopping conduc¬ 
tion, as reported in Eq. (29). We focus on a uniform- 
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FIG. 4. Temperature dependence of the electrical con¬ 
ductance for (a) N = 4-site wire and (b) N = 10-site 
wire, demonstrating tunneling behavior at low temperatures, 
ballistic-resonant activated conduction at high temperatures, 
and hopping behavior when 7 d ^ 0 , at low-intermediate 
temperatures. Parameters are eb = 0.5, v = 0.05, and 
JL = Jr = 0.2, Afi = 0.001 eV, dephasing probe condition. 



FIG. 5. Activation energy Ea, as resolved from Arrhe¬ 
nius plots as in Fig. 4, plotted against the bridge energy tB, 
demonstrating a linear relationship. Parameters are indicated 
in Fig. 4. 


symmetric system, €r = e n = 0.5—1 eV, v n ^ n ±\ = v = 
0.05 — 0.3 eV and 7 l = Jr = 0.2 — 1 eV. The tem¬ 
perature of the L and R leads, as well as the probes, is 
taken within the range T = 150 — 400 K, fixed across the 
junction. In this section we restrict simulations to the 
low voltage regime, Afi = 10 -4 — 0.01 eV. The voltage 
shifts the position of the Fermi energy at the contacts, 
Ml = Mf + eAV/2, /j,r = mf — eAV/2, but bridge ener¬ 
gies are intact. In Sec. V we study the role of a large 
applied voltage A/x ~ 0.1 — 2 eV on the transport behav¬ 
ior. The dephasing strength 7 d (probe hybridization) is 
taken uniform along the chain. We provide it in units of 
eV; when divided by h it translates to a rate constant. 

The net current flowing under an applied bias A/x is 
calculated in three steps: (i) We construct the matrices 
G and T a , then the transmission coefficient between ev¬ 
ery two terminals T a , a '(e) from Eq. (10). (ii) We include 
environmental effects e.g. with dephasing probes. This is 
achieved by calculating (numerically) the functions /„ (e) 
using the probe condition (14). (iii) We evaluate the 
net-total current from the L electrode via Eq. (7). As 


mentioned above, the probe conditions ensure that the 
current leaving the L metal is identical to the current 
arriving at the right electrode, I = II = — Ir. For de¬ 
phasing probe, moreover, 1(e) = /^(e) = —Ir{( 0- 


A. Length dependence 

The behavior of the conductance with increasing 
molecular length is displayed in Fig. 3 under a range 
of dephasing energies, 7 d = 0 — 100 meV. In the ab¬ 
sence of dephasing the conductance decays exponentially 
with molecular length for N ~ 1 — 5, in agreement with 
Eq. (27), but beyond that it is fixed (and small), inde¬ 
pendent of size. This residual component emerges due 
to incoming electrons of energy cr (occupying the tail 
of the Fermi function), which are crossing the junction 
ballistically. In contrast, at finite dephasing strengths 
and beyond N ~ 5 the tunneling behavior is overpow¬ 
ered by a hopping-ohmic contribution G^ 1 oc N. Fig. 
3 establishes two central results for hopping conduction: 
at finite dephasing strength long-range electron transfer 
follows an ohmic form, and that G oc 7^ in this region. 


B. Thermal activation 

The onset of an Arrhenius behavior at high enough 
temperatures is considered a central characteristic of hop¬ 
ping dynamics in donor-bridge-acceptor electron trans¬ 
fer processes. However, Eq. (24) reveals that in metal- 
molecule-metal geometries a thermally activated conduc¬ 
tion can take place even in the coherent regime, support¬ 
ing ballistic-resonant behavior. In fact, electrons which 
are thermally activated in the metal are more likely to 
cross the junction ballistically rather than in an ohmic 
fashion (after being scattered between probes). Indeed 
we find that a clear Arrhenius dependency in our simula¬ 
tions corresponds to ballistic conduction, while a weaker 
temperature dependency may correspond to thermalized- 
yet-off-resonant electrons, hopping across the bridge. 

Figure 4 displays this behavior. First, focus on the 
case of 7 d = 0 (red line). At very low temperatures the 
conductance does not depend on temperature at all (tun¬ 
neling contribution). At high temperatures, an Arrhenius 
behavior is observed, with an activation energy Ea, in¬ 
dicating ballistic conduction at zero dephasing. Fig. 5 
demonstrates that the activation energy Ea, the slope in 
the Arrhenius regime, logG oc —/3, linearly follows the 
bridge height cr- 

At finite dephasing strength the situation is rather in¬ 
volved. Let us first focus on the low temperature regime 
of Fig. 4. While the temperature dependence is rather 
weak, it is clear that G is increasing as G oc 7 J. This in¬ 
dicates on an activation-less dephasing assisted conduc¬ 
tance, see Eq. (29). At high temperatures and for long 
wires (see for example N = 10, T _1 = 0.03 K _1 ) the con¬ 
ductance shows an activated form, but it does not depend 
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on 7d, for y^ = 1 — 10 meV. This again points to a ballistic 
motion; the dominant role of temperature is to enhance 
the population of electrons in resonance with the bridge, 
with the ballistic component dominating hopping con¬ 
tribution. Thermally-activated hopping conduction (29) 
further shows up, but only when the dephasing strength 
is high and the bridge is short, see TV = 4, T -1 = 0.02 
K -1 and y^ = 100 meV. 

To summarize, unlike DBA-ET rates, activated con¬ 
ductance in junctions does not expose transport mech¬ 
anisms: At low temperatures, a seemingly activation¬ 
less conductance may reflect either deep-tunneling or 
dephasing-assisted conduction of deep electrons. Simi¬ 
larly, activated transport at high temperatures may cor¬ 
respond to both ballistic and ohmic components. 


C. Dephasing 



FIG. 6 . Kramers-like behavior of the electrical conductance 
with dephasing strength (dephasing probe). Parameters are 
indicated within the figure. 


The role of dephasing strength y^ on the conductance 
is explored in Fig. 6, and we observe a characteristic 
“Kramers-like” turnover behavior 82 . At low dephasing 
the conductance increases with y^ as it opens up a new 
route for electrons to cross the bridge, by hopping be¬ 
tween sites. At large dephasing strengths 7 d/tB > 2 
the conductance drops approximately as y^ 1 . In a clas¬ 
sical language one attributes this decay to the imped¬ 
iment of electrons at strong friction, resulting in an 
overdamped dynamics. This Kramers-like turnover be¬ 
havior was demonstrated in other studies of electron 
transfer 41,42 ’ 46,47 . Investigations based on Redfield 42 and 
Lindblad equations 4 ' even provided a very similar value 
for the turnover point. Related trends were also observed 
in recent explorations of environmental-assisted quantum 
energy transport in light-harvesting biomolecules 83 . 

Note that in the Kramers’ theory the barrier crossing 
rate first increases linearly with friction strength, then 
decreases as the inverse of friction. Here, we observe a 
more intricate behavior: The LBP conductance increases 
linearly with dephasing in short chains, G oc 7d, but in 
long chains the enhancement follows G oc y 2 , in the weak- 
dephasing regime. 

The quadratic dependence G oc yj is a direct result of 
the phenomenological-probe modelling of bath-assisted 
electron scatterings. Unlike a genuine electron-phonon 
scattering process, in which an electron is scattered be¬ 
tween different electronic states assisted by a phonon, 
the probes absorb an electron and emit it (possibly) at a 
different molecular state, to be re-absorbed by another 
probe. This probe-molecule-probe transfer process of 
electrons results in the yj functional form, in analogy 
with the 7 l x y# term in Eq. (27). In short chains, the 
L and R metals-to-probes scattering is significant, more 
than multiple probe-molecule-probe electron scattering, 
to construct the prefactor yy^. While the quadratic de¬ 
pendence is not necessarily physical, it is yet encouraging 
to find that the Kramers-like turnover is well reproduced 
here, in a seemingly correct value. 


D. Junction’s energetics 


To further expose the dependence of the hopping con¬ 
ductance on the junction’s energy parameters, we study 
the behavior of G with eg (set relative to the Fermi en¬ 
ergy), inter-site tunneling v, and the hybridization energy 
7 l,r- As before, we restrict parameters to the regime 
es > «,7 l , r - 

In Fig. 7 we display the conductance as a function of 
cb for TV = 8 and TV = 12 under nonzero dephasing and 
reveal a crossover from G ~ e _/3es to G ~ corre- 

e B 

sponding to the transition between a ballistic behavior 
(no length dependence) dominating at small eg, to off- 
resonance dephasing-induced hopping conduction. 

We now analyze more carefully the dependence of the 
conductance on temperature and bridge energy. We re¬ 
call that two limiting mechanisms prevail in the absence 
of dephasing effects: deep tunneling conduction and bal¬ 
listic transmission. Ballistic motion dominates in long 
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FIG. 7. Conductance as a function of the bridge energy 
tB- v = 0.05, 7 l,r = 0.2, 7 d = 10 ~ 3 (dephasing probe) and 
Ay = 0.01 eV. 
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FIG. 8 . Analysis of bath-assisted conduction, = 1 meV. (a) At low temperatures T = 150 K and for high bridges an 
activation-less transport takes place, (b) At high temperatures T = 350 K the conductance is enhanced by a thermal factor. 
Panels (al)-(a4) and (bl)-(b4) expose the ohmic-like behavior in long wires, v = 0.05, 7 l,r = 0.2, Ap = 0.01 eV, dephasing 
probe condition (full symbols), voltage probe (empty symbols). 


wires, and it carries the factor exp(— (3eb), see Eqs. (23)- 
(24), reflecting the thermal occupation of electrons in the 
metal at the bridge energy eb- Since it is non-dissipative, 
the ballistic component does not decay with distance, see 
Fig. 3. Dephasing-assisted conductance has a nontrivial 
temperature dependency due to the contribution of many 
electrons in the metal. In Fig. 8 (a) we find that at low 
temperatures the conductance is missing an Arrhenius 
activation factor, G(eb = 0.4)/G(cb = 0 . 8 ) ~ (frf ) 4 = 
16, while a different behavior takes place at high temper¬ 
atures T = 350 K, as demonstrated in panel (b). Here, in 
conjunction with an ohmic-like decay, (focus on the red 
and black lines) G{eb = 0.4 )/G{eb = 0 . 8 ) ~ 10 5 , reflect¬ 
ing an activated-thermally enhanced conductance. We 
confirm in panels (bl)-(b3) that an ohmic-like behavior 
develops with 1/G oc N. 


We also compare in Fig. 8 the operation of the de¬ 
phasing and voltage probes. We find that at this low bias 
limit (A/i = 10 meV) transport characteristics are almost 
identical under either probes: the tunneling-to-hopping 
crossover takes place at a similar molecular bridge size, 
and the values for the conductances are close. The two 
probes show some deviations only when the bridge height 
is taken relatively low eb — 0.2 eV, with the dephasing 
probe yielding lower values for conductance (larger resis¬ 
tance per site) in the hopping regime. 


E. Hopping conduction: construction of Eq. (29). 

Based on the simulations reported above, as well as ad¬ 
ditional results included in the Appendix, we suggest the 

2 2 

hopping conductance, Eq. (29), Gh ~ GoA(T)-^—jAg, 
with the dimensionless coefficient Z, positive or negative, 
see Fig 7. This expression constitutes central charac¬ 
teristics of bath-assisted conduction: Ohmic-like conduc¬ 
tion for long chains, and the enhancement of G with de¬ 
phasing strength (for weak dephasing) and with barrier 
parameters as Recall that Eq. (29) was developed 
assuming parameters in the range ")d, v, /? , A/i <C eb, 

7 < es, 

We can justify this hopping form with some qualitative 

arguments: Electrons scatter into the molecule through 

the resonance most strongly hybridized with the L metal, 

and leave from the orbital most tightly coupled to the R 

2 

lead, overall contributing the factor Ar- On the molecule, 
the probes absorb and emit electrons, dispersing them 
between molecular orbitals, thus the conductance scales 
as 7 ^, describing probe-molecule-probe transitions. The 
contribution of probe-mediated scattering processes lin¬ 
early grows with the number of molecular sites (probes), 
thus the hopping resistance grows as N. The overall tem¬ 
perature dependence in the process is rather weak; the 
hopping conductance integrates electrons thermally pop¬ 
ulating many levels in the lead, off-resonant as well as in 
resonance with molecular orbitals. 

Both dephasing and voltage probes disperse electrons 
between molecular orbitals. However, the dephasing 
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probe does so while maintaining the number of electrons 
within each energy interval fixed. The voltage probe can 
absorb/provide energy from/to electrons, and it can thus 
modify the energy distribution of electrons in the system. 
At low bias, A/z < 7 W ,/ 3 _1 , the two probes produce very 
similar conductance characteristics. 


F. Non-Markovian environment 

We implement a structured environment by using 
energy-dependent dephasing energies with the Debye- 
Drude form, Eq. (19), see Fig. 9. As expected, we 
find that the conductance in the hopping regime is re¬ 
duced due to the non-Markovianity of the bath, while 
the tunneling and ballistic regimes are unaffected. The 
effect is clearly observed at low temperatures, but it is 
significantly mitigated at higher temperatures, when the 
dominant mechanism is ballistic conduction. 



FIG. 9. Electrical conductance as a function of length using 
structured [Debye-Drude form, Eq. (19)] and unstructured 
( 7 d. is a constant) environments at high and low temperatures. 
At low temperatures the conductance is sensitive to the cutoff 
frequency O = 0.1 eV. We used v = 0.05, tB = 0.5, 7 l,r = 
0.2, A/z = 0.01 eV and 7 d = 1,10 meV as indicated, dephasing 
probe condition. For unstructured environments we used 'yd = 
1,10 meV, independent of energy. 


V. BEYOND LINEAR RESPONSE: 
CURRENT-VOLTAGE CHARACTERISTICS 


in several works 84-87 , in the weak metal-molecule cou¬ 
pling limit the applied voltage typically drops only at 
the metal-molecule contacts (symmetrically or asymmet¬ 
rically), thus the barrier maintains its equilibrium form. 
In contrast, at intermediate-strong molecule-metal cou¬ 
pling the applied bias largely drops across the molecule. 
Assuming a linear drop, the molecular barrier trans¬ 
forms from a rectangular to a trapezoidal shape, then 
to a triangular form when the bias exceeds the barrier 
height. This transition from a “direct tunneling” to 
“field emission transport” (Fowler-Nordheim tunneling) 
was explored experimentally with transition voltage spec¬ 
troscopy (TVS), when plotting log (I/V 2 ) as a function 
of l/V to extract the barrier height from the minima of 
this plot. 

Here, we assume a linear potential profile and shift the 
energies of sites n = 1, 2,..., TV accordingly, 





A/z(n — 1 ) 
TV — 1 


(31) 


Using this model, I — V characteristics at different values 
for the dephasing strength are plotted in Fig. 10, as well 
as examples for the potential profile at low, intermediate, 
and large biases, see inset. The I — V curves reveal the 
followings: Up to approximately 0.1 V a linear response 
behavior is maintained. The current is strongly enhanced 
when the voltage bias approaches a molecular resonance 
(located at ~ 0.5 eV). Beyond that, A/z >2 es, the cur¬ 
rent saturates. The effect of the thermal bath on the 
I — V characteristics is substantial. At low bias, dephas¬ 
ing effects assist electron conduction, G oc 7 ^. Further, 
in resonant situations while the current is saturated for 
7 d = 0 , it displays a negative differential conductance un¬ 
der dephasing. At large bias, the effect of dephasing on 
current is non-monotonic and deserves further attention. 

We point that a detailed TVS analysis, or the applica¬ 
tion of similar 88 or alternative 89 measures, are not par¬ 
ticularly illuminating here since we had adjusted the bar¬ 
rier, from a rectangular to a triangular shape, by hand, 
imposing the model (31). A TVS analysis or related in¬ 
vestigations are interesting for e.g., exposing the evolu¬ 
tion of the barrier height with size and contact proper¬ 
ties in different families of molecules, such as alkanethiols 
chains or aromatic thiols systems 87 . 


VI. COMPARISON TO EXPERIMENTS: ONI WIRES 


As a further application of our method, we present in 
Fig. 10 I — V traces of molecular wires of size N = 7. 
These simulations were performed using the dephasing 
probe condition since the present implementation of the 
voltage probe is applicable only in the linear response 
regime. The voltage probe technique can be extended 
beyond that, by solving sets of nonlinear equations, see 
e.g. Ref. 80 . 

Considering molecules under a large applied bias, we 
need to first determine the details of the electrostatic 
potential profile across the junction. As was discussed 


In this Section we examine experimental results us¬ 
ing the LBP method, to demonstrate the utility of 
the technique. Ref . 20 reported the resistance and 
current-voltage characteristics of conjugated oligonaph- 
thalenefluoreneimine (ONI) molecular wires of differ¬ 
ent lengths, up to 10 nm with 10 repeating units, 
measured using conducting-probe atomic force (CP- 
AFM) microscopy. Based on the length-dependent and 
temperature-dependent resistance behavior, it was ar¬ 
gued that short wires conduct via a coherent-tunneling 







12 



FIG. 10. I — V traces for a molecular bridge seven units 
long with T = 298 K, v = 0.05, and 7 l,r = 0.2 eV. The 
bridge energies shift with the applied bias in a linear manner, 
Eq. (31), see examples in the inset. Sites n = 0 and n = 8 
correspond to the chemical potentials at the contacts, rl and 
Hr, respectively. 

transport mechanism, while around 4 nm a transition 
in transport behavior, from tunneling to thermally- 
activated hopping, took place. Below, we use the data 
reported in Ref . 20 and extract relevant parameters for 
our model. We then demonstrate in Figs. 11-12 that by 
tweaking the dephasing rate and the temperature we can 
reasonably reproduce experimental results of transport 
in ONI wires, further providing some insights. 

It should be emphasized that values for resistance re¬ 
ported in Ref . 20 correspond to a monolayer with M ~ 
100; M is an estimate for the number of molecules un¬ 
der the measuring tip 20 . If we assume a linear scaling of 
the conductance, G(M) ~ MG( 1), the resistance should 
obey R(M) ~ M~ 1 R{1). Certainly, this is an approxi¬ 
mation which is not always justified; the conductance of 
a single-isolated wire G(l) may significantly deviate from 
the monolayer conductance scaled-down, G(M)/M, due 
to cooperative effects, inter-wire couplings and substrate- 
mediated coupling. These effects may increase or de¬ 
crease the conductance per wire 90-94 . A careful study of 
such effects with the LBP method is of interest. Here 
instead we use a naive-linear scaling approximation: We 
take resistance values reported in Figs. 6 and 8 of Ref . 20 
and multiply them by M = 100; see the symbols (+) in 
Figs. 11 and 12. Below we explain which parameters in 
our modelling are affected by this simplified scaling, and 
if so, the impact on analysis. 

To uncover model parameters from Ref. 20 , cb, v, Jl,r, 
and 7 d, we first determine the unit length of the ONI 
wires to be a ~ 1 nm. We then obtain our parameters 
by following three steps: 

(a) Bridge energetics. In the tunneling regime N = 
1 — 3, G ~ e~ KNa , with k = 2.5 nm -120 . Using the re¬ 
lation k = | In 1 7 ^ | [see Eq. (27)], we resolve the ratio 


£b/v = 3.4. To obtain the bridge height cb and inter-site 
energies, separately, we note that the activation energy 
in ONI wires was approximated by Ea ~ 0.55 eV from an 
Arrhenius plot 20 . Furthermore, in Fig. 4 we found that 
while the activation energy linearly follows the bridge 
height, the slope deviates from unity, cb ~ IAEa for 
7 d ~ 5 meV. We thus select the value cb = 0.8 eV, and 
immediately receive v = 0.22 eV. Note that these esti¬ 
mates are not affected by the scaling employed. Note also 
that the ballistic conductance is very sensitive to temper¬ 
ature near room temperature, and therefore is also very 
sensitive to the bridge height/activation energy. Even 
when using our “corrected” cb = 0.8 value for bridge 
energy, we still observe excess ballistic conductance near 
room temperature, and can thus only match experimen¬ 
tal results at a lower, “effective temperature”. 

(b) Hybridization energy. We resolve Jl,r from Eqs. 
(26)-(27) (using Go = e 2 /itH to accommodate spin de¬ 
gree of freedom). We assume that this contact energy 
is identical at the two interfaces, and employ = 
s/GjG u x e K,Na . From the measured exponent k = 2.5 
nm -1 we receive 7 „ ~ 0.018 eV. This estimate will vary 
if we adjust the scaling for G(M), see discussion below. 

(c) Dephasing strength. We approach this parame¬ 
ter, characterizing electron-vibration interaction in the 
molecule, by testing different values (uniform along the 
wire) in the range 7 d = 0 — 100 meV. In Fig. 11 we 
show that very small dephasing values 7 ^ = 0 — 0.2 meV 
reasonably reproduce both the exponent in the tunnel¬ 
ing regime and the slope in the hopping-ohmic domain, 
assuming the junction’s temperature is in the range 280- 
300 K. As mentioned above, a very good agreement is 
reached by slightly reducing the temperature; alterna¬ 
tively, higher bridge heights could be tested. 

With these molecular parameters, we turn now to the 
Arrhenius plot of ONI wires, as reported in Ref. 20 . Fig¬ 
ure 12 depicts (+) the resistance as a function of inverse 
temperature for an IV = 7 ONI wire, taken from Ref. 20 , 
again linearly scaled down to the single molecule value. 
We find that our results agree fairly well with experimen¬ 
tal data under a coherent transport mechanism, and we 
resolve the activation energy Ea ~ 0.4 eV from simu¬ 
lations, compared to the experimental value Ea ~ 0.5 
eV. We try to reproduce this data using the parame¬ 
ters resolved in the discussion above, but find out that 
while at room temperatures and above, the experimen¬ 
tal data is reasonably captured by our simulations with 
7 d ~ 3 — 4 meV, at lower temperatures the experiment 
is reproduced with 7 ^ = 0 meV. A possible explanation 
for this behavior is that the dephasing rate in fact de¬ 
pends on temperature, 7 d(T), and its value should be 
taken small 7 ^ -C 1 meV at low temperatures, T = 240 
K, while at room temperature it reaches values identified 
above. 

Our analysis is based on the linear scaling assumption. 
It was demonstrated in Ref . 90 that this holds only beyond 
a certain molecular-island size of a few tens of molecules. 
Scaling the experimental results of Ref . 20 by M / 100 
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would alter our estimate for and 7 recall that G oc 
in the tunneling regime and G oc 7 J in the hopping 
limit, see Sec. IV. As a limiting case, if reported values 
in Ref . 20 were to correspond to a single molecule rather 
than to M ~ 100 molecules, our simulations would fit 
the data with the same values for the bridge energetics, 
eb = 0.8 eV, v = 0.22, but with (factor of 10 larger) 
7 l = 1r = 0.18 eV and 7 ^ ~ 2 meV. These values should 
serve as an upper estimate for dephasing strength. 

Concluding our observations, the probe technique un¬ 
der weak dephasing as a tweaking parameter can capture 
the turnover behavior between superexchange and hop¬ 
ping regimes in molecular wires, given the correct system 
parameters as inputs. We also see from Fig. Il(a2), as 
well as in Fig. 12, that in order to match the experiment 
across all temperatures and lengths, the tuning parame¬ 
ter 7 ^ must depend on temperature, and that one should 
take into account corrections to the bridge height as it 
relates to the activation energy. 


+ ONI data 

Empty squares: simulations 
N = 7, 7 ^ = 0, 0.2, 0.4 meV ^' 


N = 3 

o .O.O.O.O.O . O.O.G | 

3 32 3 A 3i> 343 

1000/T [if 1 ] 

FIG. 12. Analysis of experimental results for the temperature 
dependence of conductance in ONI wires. Experimental data 
(_l_)2° f or tv = 7 j s compared to simulations with 7d=0, 0.2, 
0.4 meV (empty squares, top to bottom). For reference, we 
also show the N = 3 case (o) in which the resistance is inde¬ 
pendent of dephasing in the range 7d ~ 0 — 5 meV (consistent 
with the experiment). Other parameters are the same as in 
Fig. 11. 
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FIG. 11. Analysis of experimental results. We extracted data 
for the resistance of an ONI monolayer with M ~ 100 from 
Fig. 6 of Ref. 20 , and multiplied it by M = 100 to approximate 
resistance per molecule. Top: a single ONI wire (length de¬ 
noted by P here, rather than N, to eliminate confusion), (a)- 
(b) Experimental data (+) 20 as a function of molecular length, 
compared to probe-method simulations with 7^ = 0, 0.2, 0.5 
meV, top to bottom (dotted lines with empty squares) assum¬ 
ing (al)-(a2) T = 300 K and (bl)-(b2) T = 280 K. Panels 
(al)-(bl) illustrate the tunneling region for IV = 1 — 4, panels 
(a2)-(b2) uncover the ohmic behavior for N > 5. Dephasing- 
probe simulations were performed with eb = 0.8 eV, v = 0.22 
eV, 7_L,i?=0.018 eV, applied voltage Ap = 0.01 eV and de¬ 
phasing strengths in the range 7d = 0 — 1 meV. 


VII. CONCLUSIONS 


We demonstrated that the Landauer-Biittiker probe 
technique can be used to simulate electronic conduction 
in molecular junctions under the influence of a ther¬ 
mal environment. The LBP approach can be incorpo¬ 
rated with little effort into the commonly used Landauer- 
NEGF approach, to phenomenologically introduce bath- 
assisted site-to-site hopping conduction. 

Our simulations revealed different transport mecha¬ 
nisms: In the absence of dephasing effects we identi¬ 
fied deep-tunneling conduction in short wires and ballis¬ 
tic transport at high temperatures, a result of resonant 
transmission. Under finite dephasing strengths, hopping 
conduction was the dominant transport mechanism in 
wires of N > 4, showing temperature dependence which 
was generally weaker than the Arrhenius activation fac¬ 
tor in the explored range of parameters. 

Principal observations are: (i) Length-dependence in¬ 
vestigations of molecular conduction immediately pin¬ 
point the tunneling-to ohmic turnover. In contrast, ac¬ 
tivated conduction (or more generally, temperature de¬ 
pendent conduction) may show up in both the coherent 
(ballistic) and incoherent-hopping regimes, (ii) Kramers- 
like turnover behavior manifests itself under the LBP 
method. Weak dephasing effects promote hopping con¬ 
duction, but dephasing becomes detrimental to transport 
when 7 d > eb- (hi) Structured environments can be 
implemented within the LBP method, by using energy- 
dependent dephasing coefficients, (iv) At large voltage 
biases, dephasing effects realize a negative differential 
conductance; this behavior is missing in the coherent 
limit, (v) A closed-form expression for the hopping con¬ 
duction was constructed, Eq. (29), within stated parame¬ 
ters. (vi) The LBP method provided a semi-quantitative 
match for length dependent resistance in ONI wires re- 
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producing the tunneling to hopping crossover, as well as 
the resistance’s temperature dependence- as long as we 
allow the dephasing strength to become temperature de¬ 
pendent. 

Biittiker’s probes offer an easily implemented mean to 
introduce dephasing and inelastic effects into molecular 
electronic applications, yet, obviously they do not abso¬ 
lutely correctly emulate physical process e.g., electron- 
vibration interaction. While overall the LBP approach 
ensures the conservation of current across the molecule, 
even within a small energy interval (dephasing probe), 
one should remember that the probes are metal elec¬ 
trodes, collection of Fermi sea electrons which themselves 
participate in the transport process. Essentially, the 
quadratic dependence of the hopping conductance on the 
probe hybridization ( Gh oc 7 2 ) reflects probe-molecule- 
probe electron scatterings, with an incoming electron 
from probe n arriving at the molecule, leaving it to 
probe n'\ complementary processes ensure conservation 
of charge and energy, as necessary. In contrast, micro¬ 
scopic modelling of electron-vibration interactions should 
yield to lowest order a hopping conductance linear in the 
dephasing rate, possibly modifying Eq. (29), 7 ^ —> "fd/P- 
Thus, while the probes reproduce the hopping-ohmic na¬ 
ture of bath-assisted conduction, the predicted dephasing 
and temperature dependence should be taken with cau¬ 
tion. To understand the correspondence of LBP results 
with physical-microscopic modelling, one should compare 
our simulations to projection-operator approaches and 
first-principle techniques 36 . 

In future studies we will examine the suitability of the 
LBP technique to describe the phenomenology of elec¬ 
tronic function under environmental effects, e.g. consider 
molecular thermoelectricity under thermal effects 52,95 
and the development of a diode behavior in symmetric 
molecules, induced by the surroundings 96 . 
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APPENDIX A: SUPPORTING SIMULATIONS FOR 
HOPPING CONDUCTANCE. 

In the main text we established Eq. (29) by explor¬ 
ing the dependence of the hopping conductance on cb, 
jd, temperature and length. Here we complement this 
and examine the role of inter-site tunneling v and metal- 
molecule hybridization on Gh- In Fig. 13 we confirm 
that in the coherent-tunneling regime G oc v 2N ~ 2 , but at 


finite dephasing and for large N we obtain G oc v 2 for 
independent of molecular length (not shown). 


.-N = 4,T = 200R,7 d = 0 
- » -TV = 12, T = 300R,7rf = 10 meV 



FIG. 13. Electrical conductance vs. inter-site tunneling v. In 

1 / \ 2JV 

the coherent tunneling limit G oc ^ I ^ ) .In the hopping 

regime G oc v 2 . Parameters used are es = 0.5, v = 0.05, 
7 „ = 0.2, Afj, = 0.01 eV. 



FIG. 14. Electrical conductance as a function of the hy¬ 
bridization energy, 7 l = 7 r- When 7^ = 0 we resolve tun¬ 
neling conductance for IV = 3, G oc ^LfR (A) and ballistic 
motion for N = 7 when 7 „ cb (o). The conductance depen¬ 
dency on 7 „ diminishes upon increasing dephasing strength. 
We used T = 300 K, es = 0.5, v = 0.05, Afj, = 0.01 eV, 
dephasing probe simulations. 

The hybridization energy 7 ^ is expected to only mildly 
influence hopping conduction, as the resistance should 
be determined by the wire itself - electron scattering be¬ 
tween probes- rather than by the interface of the molecule 
with the contacts. This behavior is confirmed in Fig. 14: 
Tunneling conductance ( N = 3, 7 ^ = 0) grows as 7 l 7 _r, 
see Eq. (27), but in long chains with 7 ^ yf 0 the effect of 
the hybridization energy on G is insignificant. 
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